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ABSTRACT 

We present WHT and VLT spectroscopy of MM Ser, the optical counterpart to Ser X- 
1. We deblend the red spectra of the two close stars identified by Wachter (1997) and 
show that the brighter of the two is responsible for the Ha and He 1 emission, hence 
confirming that this is the true counterpart of the X-ray source. We also identify several 
He 11 and N in lines in the blue spectrum. The isolated emission lines are all remarkably 
narrow, with FWHM 200-300 km s _1 . The Bowen blend has structure suggesting that 
the individual components are also narrow. These narrow lines could be from the disc 
if the binary inclination is quite low, or they could come from a more localised region 
such as the heated face of the companion star. Several interstellar lines are detected 
and indicate that the reddening is moderate, and consistent with the neutral hydrogen 
column density inferred in X-rays. 
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1 INTRODUCTION 

The low-mass X-ray binary (LMXB) Ser X-l has been 
known as an X-ray source from the early days of X-ray as- 
tronomy. It was discovered in 1965 (Friedman, Byrom & 
Chubb 1967). Its optical counterpart has proven more elu- 
sive, however. When an accurate (1 arcmin) position for the 
X-ray source became known (Doxsey 1975), Davidsen (1975) 
suggested that the optical counterpart was an ultraviolet 
excess object with B ~ 18.5. Subsequently Thorstensen, 
Charles & Bowyer (1980) used images obtained in better 
seeing conditions to show that this counterpart was actually 
two coincident stars (DN and DS) separated by 2.1arcsec. 
The southern one, designated MM Ser, was by far the 
brighter in the ultraviolet. The detection of an optical burst 
from DS, simultaneously with an X-ray burst (Hackwell et 
al. 1979), confirmed that the X-ray source was associated 
with DS. More recently, however, Wachter (1997) has shown 
that DS is itself two unresolved stars (DSe and DSw), sepa- 
rated by only 1 arcsec. Wachter suggested that the brighter 
of the two stars, which is bluer, might be the true optical 
counterpart. Silber (1998) has analysed rapid photometry 
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which supports this, indicating that the brighter of the two 
stars is significantly variable. No convincing periodicity was 
found in these data, however. 

Spectroscopy of MM Ser was obtained at several epochs, 
before it was resolved into two stars. Thorstensen et al. 
(1980) obtained blue spectroscopy of both DS and DN. Both 
showed H/3 absorption, but DS in addition had a He 11 4686 A 
emission line. Cowley, Hutchings & Crampton (1988) con- 
firmed the presence of narrow He 11 emission. Their spectrum 
shows no H/3 absorption, but possible N ill 4640 A emission. 
Shahbaz et al. (1996) reported an almost featureless contin- 
uum spanning the full optical bandpass, albeit at a very low 
resolution which would be insensitive to weak emission fea- 
tures. An absorption feature around 5900 A was attributed 
to a G star secondary. 

We report on further spectroscopy of MM Ser. Our pri- 
mary goal was to separate the spectra of DSe and DSw and 
hence confirm which is the true optical counterpart to Ser X- 
1. Since the most recent published spectroscopy dates from 
1988 we also were able to obtain a higher quality spectrum 
using modern instrumentation and hence study the proper- 
ties of MM Ser in more detail than previously possible. 
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2 OBSERVATIONS 
2.1 WHT spectra 

Spectra of MM Ser were obtained from the Observatorio del 
Roque de los Muchachos on La Palma on the nights of 2001 
July 8 and 10 using the 4.2 m William Herschel Telescope 
(WHT) with the ISIS dual-arm spectrograph. The seeing 
was typically 1.0-1.5 arcsec on July 8 and ~ 0.8arcsec on 
July 10. 

The slit width was set to 0.83 arcsec. The slit was 
aligned to pass through the two blended stars (DSe and 
DSw) identified by Wachter (1997). The blue arm was used 
with the EE VI 2 CCD and the R600B grating to yield a 
resolution of 1.9 A. Note that the blue spectrum is heav- 
ily vignetted at either end; the useful coverage spans 3720- 
5260 A. On July 8 the red arm was used with the TEK4 
CCD and the R600R grating for a resolution of 1.5 A. On 
July 10 the R316R grating was instead used giving a reso- 
lution of 3.1 A. A log of the spectra obtained of MM Ser is 
given in Tabled 

Basic image processing (bias correction and flat fielding) 
was done using standard iraf 1 techniques. Unfortunately no 
suitable red arm flat field was obtained on July 10, so it was 
necessary to use the July 8 one which was obtained with a 
different grating. 

Initially a combined spectrum of both stars was ex- 
tracted using the iraf implementation of optimal extrac- 
tion (Home 1986; Marsh 1989). As there are actually two 
blended stars at the site of MM Ser, this will not truly be 
an optimal extraction in this case. The spectra extracted in 
this way were still far superior to a non-optimal extraction, 
however. 

Wavelength calibration was interpolated relative to ex- 
posures of a CuNe+CuAr lamp and checked against night 
sky lines. The rms scatter among the latter was 0.25 A in 
the blue and 0.14 A in the red, with no systematic offset in 
either case. 

Flux calibration was applied relative to the spectropho- 
tometric standard BD+33°2642 (Oke 1990); this was also 
used to correct for Telluric absorption features. As an 
0.8 arcsec slit was used significant uncertainty will be in- 
troduced in the flux calibration due to slit losses, especially 
as the slit was not aligned at the parallactic angle due to 
the need to follow the line of centres of the blended star. 
Consequently, the flux calibration of the spectra is only ap- 
proximate. 

The average extracted spectra (a composite of the two 
blended stars) is shown in Fig. The presence of strong 
emission lines in the combined spectrum suggests that the 
brighter, bluer star is the optical counterpart to MM Ser as 
suggested by Wachter (1997). To confirm this we also used 
the spectral deblending method of Hynes (2002) to extract 
spectra of the two stars separately. 



IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 




3800 4000 4200 4400 4600 4800 5000 5200 
Wavelength (A) 



2.0 



E 1.5 



1 .0 



0.5 



(b) 



Hel 



Ha 




DIB 




DIB 



NaD 



6000 6200 6400 6600 6800 7000 7200 
Wavelength (A) 

Figure 1. a) Combined blue spectrum of the two stars, binned 
to 3.5 A pixels. This will be dominated by MM Ser with the other 
star contributing less than 20 percent of the light, b) Combined 
red spectrum. Identified emission and interstellar absorption fea- 
tures are marked. Telluric correction has been applied, but re- 
gions which may be contaminated by residuals from this process 
are marked ffi. The Telluric feature around 6300 A is much weaker 
than the longer wavelength features, so this region is probably not 
significantly contaminated. 



2.2 VLT Spectra 

Three further spectra were obtained with the FORS2 spec- 
trograph at the Very Large Telescope (VLT) on 2002 May 
18. The 1400V grism was used with an MIT/Lincoln Labs 
mosaic CCD and an 0.7 arcsec slit. The seeing was around 
1.2 arcsec. Further details are given in Table 

Images were debiased and fiatfielded using Starlink 
CCDPACK routines. As no suitable PSF template was avail- 
able, and the contribution of the fainter star appeared very 
small, we only performed a single object 'optimal' extraction 
with iraf as above. The images were marred by a number of 
bright spots, typically up to 20-30 pixels across. These orig- 
inate from enhanced dark current in regions that had been 
heavily saturated during earlier observations. We examined 
the saturated images to identify which regions of the spec- 
trum might be affected. Fortunately, there were only a few 
mildly saturated regions coincident with the location of our 
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Table 1. Log of spectroscopic observations of MM Ser. 
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Figure 2. Average, normalised VLT spectrum of MM Ser. Iden- 
tified features are marked. Unidentified features appear to be 
present at 4605 and 4945 A. These are repeatable, and not ob- 
viously artifacts. 



Figure 3. Average WHT red spatial profile showing MM Ser 
and the blended star. The solid line shows the data, oversampled 
and averaged over wavelength after removing the curvature of the 
spectrum. The dashed lines show the corresponding averages of 
the two modelled profiles. 



spectra, and none of these correspond to apparent features in 
the extracted one- dimensional spectra. We therefore believe 
that this problem does not compromise our dataset. Wave- 
length calibration was applied with respect to an arc image 
taken during the day. The wavelength was checked with re- 
spect to sky lines (mainly [Ol] 5577 A). There was an offset 
of +0.21 A (13kms _1 at 4686 A) which was corrected. This 
is within the expected flexure for FORS2 (Szeifert 2002). 
Given the large slit losses due to poor seeing, no flux cali- 
bration of this spectrum was attempted. The final average 
spectrum is shown in Fig. |2] 



3 IDENTIFICATION OF THE COUNTERPART 

To separate the two stars we applied the deblending algo- 
rithm of Hynes (2002) to the July 10 red WHT spectra. 
These had the best spatial resolution and were well ex- 
posed, and were the only spectra suitable for this work. 
The stars were not fully resolved (Fig. [3J, but were suf- 
ficiently separated that they could be identified and de- 
blended. We defined the spatial profile with a base Voigt 
function, with width allowed to vary with wavelength, to- 
gether with a wavelength independent numerical correction. 
The deblended spectra are shown in Fig. 0] 

It is clear that the emission features detected do orig- 
inate from the bright star. Ha is particularly striking, but 
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Figure 4. Average of the two deblended red spectra of MM Ser 
and the nearby star from 2001 July 10. It is clear that the emission 
lines, particularly Ho, are associated with the brighter of the two 
stars. 
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He I 6678 A is also recognisable. These features are clear in 
both individual spectra as well as in the average. This pro- 
vides further confirmation of earlier suggestions (Wachter 
1997; Silber 1998) that the brighter of the two stars, DSe in 
the terminology of Wachter (1997), is the LMXB. 

Comparing the two components, at the epoch observed 
the fainter star contributes only about 15-20 percent of the 
total light in the 6000-7000 A band. For comparison, the 
magnitudes given by Wachter (1997) correspond to contri- 
butions of 19 percent at V and 25 percent at R. 



4 THE OPTICAL SPECTRUM OF MM Ser 
4.1 Emission lines 

The only emission lines clearly detected in the WHT blue 
spectrum are He II 4686 A and the N ill Bowen blend around 
4640 A; there are no strong Balmer emission lines in the blue. 
The Bowen blend contains at most weak C ill emission. Mov- 
ing to the red, the WHT spectrum is dominated by Ha, with 
He I 5875 and 6678 A also present. The VLT spectrum par- 
tially overlaps the WHT blue-side, but also fills in coverage 
around 5500 A. It shows the same blue lines, of N ill and He II 
4686 A, with a similar relative strength. The He II 5411 A line 
is also clearly present. Weak H/3 emission does appear to be 
present, on the blue-side of a broad absorption trough. This 
is a ubiquitous feature in both black hole and neutron star 
LMXB spectra, e.g. Buxton & Vennes (2003); see also Soria 
et al. (2000) and references therein for other black hole ex- 
amples and Hynes et al. (2001) and Casares et al. (2003) for 
other neutron star cases. If it is indeed associated with H/3, 
the large redshift (~ 1000 km s -1 ) remains unexplained. It is 
possible that the apparent redshift is only due to filling in of 
the blue absorption by emission, however. Such absorption 
troughs are usually attributed to the inner optically thick 
accretion disc (e.g. Soria et al. 2000), but a large redshift 
would not then be expected. Apparent features are seen at 
4605 and 4945 A. They do not appear to be artifacts, but 
there is no obvious identification and no strong features are 
seen in other LMXBs at these wavelengths. Given the in- 
ferior quality of the WHT spectra, they would probably be 
undetectable there. Properties of lines measured in all of the 
spectra are summarised in Table |2] 

All of the lines are remarkably narrow, as can clearly 
be seen in Fig. |S] where the line profiles in velocity space 
are collated. The FWHM listed in Table H are based on 
single Gaussian fits to the unbinned spectra. The contribu- 
tion from the instrumental resolution has been corrected for, 
but is never dominant. For comparison the FWHM of Hell 
in 4U 1822-371 (Casares et al. 2003) is FWHM 750 kms" 1 , 
and the lines in 4U 1735-444 and 4U 1636-536 respectively 
have FWHMs of 660 kms" 1 and 1220 kms" 1 (Casares et al. 
in preparation). Such narrow emission lines are thus rather 
atypical for an LMXB. If they are disc lines, then the nar- 
rowness could indicate either a long period, and hence a 
large disc with a low outer velocity, or a low inclination, 
or a combination of both. Without fully resolving the pro- 
files and confirming that they are disc lines, it is of limited 
value to model them. Applying some simple models to the 
Hell 4686 A line, however, suggests that for typical neutron 
star LMXB parameters, if the period is less than 1 day then 
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Figure 5. Line profiles of the main, identified, non-blended lines. 
All are rather narrow and appear single-peaked, at least to the 
limits of the resolution. Spectra have been binned to approxi- 
mately one pixel per resolution element in each case. They have 
been continuum subtracted and normalised to the same total line 
flux, so that the shapes of the profiles can be compared. 



i < 20° . A longer period allows a somewhat higher inclina- 
tion; for 10 days then i < 45° . These are not very robust esti- 
mates, however, as when system parameters are known, ap- 
parent sub-Keplerian velocities are often inferred from disc 
line profiles (Marsh 1998). Alternatively the lines may not 
be dominated by the disc; they could come from the com- 
panion star or stream-impact point, or some other region 
(as in XTE J2123-058; Hynes et al. 2001), in which case 
the narrowness would be more natural. In this case, unless 
the inclination were low, we would expect the lines to move. 
There is clearly no large velocity difference between the two 
epochs, and the Hell 4686 A velocities measured are consis- 
tent within errors, especially allowing for some uncertainty 
in the absolute velocities. The WHT spectra are individu- 
ally too poor to search for shorter timescale radial velocity 
variations. The three VLT spectra are better suited. They 
are equally spaced and the 5577 A night sky line is stable to 
±0.01 A between them. We measure Gaussian fit velocities 
for Hell 4686 A in the individual spectra of 105±13kms _1 , 
1 11 ± lOkms -1 , and 133±9kms~ respectively. These sug- 
gest a small drift (over about an hour), but this is not signifi- 
cant at the 90 percent confidence level and more observations 
would be needed to confirm any motion. 
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Table 2. Measured line properties. 



Date Line LSR Velocity EW FWHM 

(kms- 1 ) (A) (kms" 1 ) 

2001 July 8-10 N III 4634,41,42 1.9 ± 0.4 
(WHT/ISIS) Hen 4686 A +127 ±11 2.1 ±0.2 200 ± 30 

Hi 6562 A (Ho) +90 ± 6 2.1 ±0.1 210 ± 20 

He I 6678 A +165 + 40 0.5 ±0.1 280 ± 80 

2002 May 18 N III 4634,41,42 1.6 ± 0.4 
(VLT/FORS2) Hen 4686 A +112 + 6 1.4 + 0.1 240 + 20 

Hen 5411 A +42 + 16 0.6 + 0.1 320 + 40 
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Figure 6. The spectral region containing the Nm 4634/41/42 A 
blend and He II 4686 A. All annotations indicate the laboratory 
wavelength of expected lines. Hen is clearly somewhat redshifted, 
and Nm may be as well. 



The region around the N ill lines and He II 4686 A is 
expanded and compared in Fig. |S| The N III blend appears 
to show repeating structure with a strong peak correspond- 
ing to the 4641/42 A lines and a weaker peak for 4634 A. 
There is no evidence for Cm 4647/50/51 emission. That 
these components are resolved suggests that the individual 
components are narrow, like the isolated lines. The ratio of 
the Bowen and Hell equivalent widths is ~ 0.8, typical of 
Galactic plane LMXBs which have values of 0.5-1.0 (Motch 
& Pakull 1989). 

4.2 Interstellar features 

Our spectra show a number of interstellar features; NaD 
lines, diffuse interstellar bands (DIBs) and perhaps Ca H and 
K lines. We can use these to constrain the amount of red- 
dening toward MM Ser, although since these measures de- 
pend on measurements of equivalent widths, it is only valid 
to use measures from the deblended red spectrum. Shahbaz 
et al. (1996) suggested that the strong NaD feature in the 
spectrum of MM Ser indicated a G-type companion star. 
Since the spectrum shows strong diffuse interstellar bands 
(DIBs), we instead feel that it is more likely that the NaD 
absorption is also of interstellar origin. Call H and K lines 
(3933,3968 A) are also marginally detected in the blue spec- 



trum, which should be dominated by the accretion light; 
hence these favour an interstellar origin. The interstellar ab- 
sorption is therefore likely large. 



The NaD doublet can be used as a reddening indicator, 
but is only reliable for small reddenings (E(B — V) <0.5), 
as it is prone to be saturated at larger values (Munari & 
Zwitter 1997). We measure an EW for the NaD 5890 A line 
of 1.6 ± 0.1 A from the deblended spectrum of DSe. This is 
larger than any values given by Munari & Zwitter (1997), 
which saturate at a maximum value of about 1.2 A. There 
is a spread in the saturation values, however, as it depends 
on the relative velocities of the clouds contributing to the 
absorption. A larger maximum NaD EW is seen when mul- 
tiple resolved components are present as this partly circum- 
vents saturation. Since all of the stars considered by Munari 
& Zwitter (1997) are at < 3kpc and MM Ser has an esti- 
mated distance of 8.4 kpc (Christian & Swank 1997), it is 
likely that a larger EW is possible in the latter, since NaD 
absorption will be spread over a wider range of velocities 
and hence be less affected by saturation. We can still use 
the NaD strength as a lower limit for the reddening (if it 
is of interstellar origin), since the asymptotic linear relation 
between EW and E(B — V) at low reddenings defines the 
limiting case of unsaturated absorption. The relation given 
by Munari & Zwitter (1997) will then define a lower limit of 
of E(B — V) > 0.4. We note that the ratio of the lines is also 
an indicator of reddening (Munari & Zwitter 1997), since 
the Na D 5896 A line is intrinsically weaker and so saturates 
more slowly. At low optical depths the ratio is 2.0, but this 
decreases to 1.1 as the reddening increases. We measure a 
ratio of 1.3 ± 0.1 suggesting that at least some components 
are saturated and hence that the reddening is high, signifi- 
cantly larger than the lower limit. 



A relatively high reddening is supported by the strong 
DIBs; from the DIB at 6203 A we measure an EW of 
270+ 130 mA, implying E(B-V) = 0.8 + 0.4 (Herbig 1975). 
The DIB at 6283 A may also be strong, but will not be reli- 
able due to Telluric contamination. This value is consistent 
with the column density of Ah = 0.5 x 10 21 cm -2 estimated 
by Christian & Swank (1997), since the latter would imply 
E(B - V) = 0.6 - 1.1 (e.g. Bohlin et al. 1978; Predehl & 
Schmidtt 1995). We conclude that MM Ser is moderately 
reddened, with all indicators consistent with a reddening 
E(B-V) ~ 0.6 - 1.1. 
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5 CONCLUSIONS 

We have performed spectroscopy of MM Ser, the optical 
counterpart to Ser X-l. We resolve the two close components 
identified by Wachter (1997) and confirm earlier suggestions 
that the brighter of the two is an emission line source and 
hence the true counterpart. The spectra reveal emission lines 
of Hi, He i, Hen, and Niii, as well as interstellar features. 
The emission lines are unusually narrow. Since all lines show 
similar widths, the most likely explanation seems that they 
are disc lines, but that the binary inclination is very low. Al- 
ternatively they could also originate from a more localised 
region such as the companion star or stream-impact point, 
but we would then require that the disc not dominate any of 
the optical lines. Further time-resolved observations would 
be needed to search for motion of the lines and discriminate 
between these possibilities. Higher resolution spectroscopy 
would also be beneficial to resolve the narrow line profiles 
better and study the substructure of the N ill blend. Based 
on the strong interstellar lines, we have estimated a redden- 
ing consistent with estimates based on the X-ray derived 
neutral hydrogen column density of E(B — V) = 0.6 — 1.1. 
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